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19 .  ABSTRACT 

A  measurement  of  the  acoustic  ambient  arriving  from  a  horizontal  direc¬ 
tion  along  with  total  acoustic  intensity  spectra  allows  one  to  infer  both  the 
total  directional  spectra  and  some  physical  characteristics  of  the  sources  of 
"sea  surface  sound."  A  long-term  measurement  of  these  two  quantities  was  made 
at  high  frequency,  i.e.,  8  kHz  to  64  kHz,  in  the  Tongue  of  the  Ocean,  The 
Bahamas.  The  horizontally  directed  ambient  was  measured  using  vertically 
oriented  line  arrays  and  was  observed  for  wind  speeds  ranging  from  1  to  30 
knots,  The  resulting  data  base  was  used  to  estimate  the  statistics  of  aniso¬ 
tropic  "noise  gain"  relative  to  the  isotropic  "noise  gain."  Differences  in 
the  functional  dependence  and  residual  statistics  were  found  for  two  cases: 
whitecaps  present  and  not  present.  The  relation  of  these  results  to  the  total 
directional  spectra  and  a  model  of  the  near- surface  distribution  of  acoustic 
sources  are  discussed. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


SUMMARY 


This  report  describes  the  results  of  a  data  analysis  task  that  estimated 
the  acoustic  ambient  directional  spectra  caused  by  sea  surface  sources  of 
sound  from  measurements  of  the  horizontally  directed  acoustic  ambient.  A  the¬ 
oretical  discussion  shows  that  the  near-horizontal  value  of  the  directional 
spectra  is  a  sensitive  measure  of  the  entire  directional  spectra  resulting 
from  a  diffuse  source  field  consisting  of  surface  dipoles  and  a  subsurface 
layer  of  monopoles.  An  extended  measurement  was  undertaken  of  the  horizontal 
value  of  the  array  gain  of  three  vertical  arrays  covering  the  frequency  range 
of  8  to  32  kHz.  Measurements  of  the  array  gain  and  the  total  acoustic  inten¬ 
sity  were  made  over  a  1-year  period  along  with  a  nearby  measurement  of  the 
wind  speed  at  a  10 -m  height.  These  measurements  took  place  in  The  Tongue  of 
the  Ocean  in  The  Bahamas,  at  a  depth  of  125  m  in  1830  m  of  water.  The  results 
of  the  analysis  indicated  a  clear  wind  speed  dependence  in  the  horizontally 
steered  array  gain  resulting  from  a  wind  speed  dependence  in  the  acoustic 
ambient  vertical  directional  spectra.  The  measured  acoustic  sound  pressure 
level  spectra  were  shown  to  agree  well  with  previous  measurements.  The 
observed  array  gain  was  found  to  dichotomize  into  categories  of  either  white- 
caps  present  or  whitecaps  not  present,  with  distinctly  different  behavior  in 
each  category.  With  whitecaps  present  the  array  gain  was  found  to  be  greater 
than  the  directivity  index  of  the  array  because  the  ambient  pressure  field  was 
dominated  by  energy  arriving  nearly  vertical  to  the  array.  This  is  character¬ 
istic  of  a  source  field  of  surface  dipoles,  presumably  caused  by  whitecap  gen¬ 
eration  of  both  spray  and  an  acoustically  thin  layer  of  bubbles.  Surface 
dipoles  cause  a  maximum  array  gain  condition  which  was  observed  in  the  data  by 
examining  the  statistical  distribution  of  array  gain  values.  This  distribu¬ 
tion  was  found  to  be  skewed  to  lower  values  of  gain,  indicating  the  existence 
of  a  physical  maximum  in  array  gain  caused  by  dipole  sources.  When  whitecaps 
were  not  present  the  array  gain,  relative  to  the  directivity  index,  was  lower 
than  in  the  case  with  whitecaps  present.  While  there  was  significant  variab¬ 
ility  in  the  array  gain  fer  a  given  wind  speed,  the  array  gain  in  general  fol¬ 
lowed  a  power  law  model  with  a  1.9  exponent.  The  reduction  in  array  gain  was 
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caused  by  an  increase  in  Che  horizontal  values  of  the  vertical  directional 
spectra  relative  to  the  vertical  value.  It  is  shown  that  this  situation  could 
be  caused  by  an  acoustically  thick  sublayer  of  monopoles,  presumably  bubbles. 
If  this  model  is  correct,  the  depth  scale  of  the  vertical  distribution  of  the 
bubble  density  is  inversely  proportional  to  wind  speed.  However,  it  is  clear 
from  the  variability  of  the  array  gain  at  a  given  wind  speed  that  wind  is  only 
one  of  several  causes  of  the  change  In  the  sublayer  depth  scale.  The  data 
show,  at  least  for  this  high  frequency  range,  the  wind  speed  dependence  of  the 
acoustic  ambient  near-horizontal  directional  spectra  inferred  from  the  array 
gain  measurements.  It  is  shown  that  total  acoustic  sound  pressure  level  spec¬ 
tra  is  a  better  predictor  of  array  gain  than  is  the  wind  speed.  This  last 
result  indicates  that  when  whitecaps  are  not  present  the  wind  speed  is  not  a 
complete  indication  of  the  acoustic  state  of  the  sea  and,  furthermore,  that 
total  acoustic  sound  pressure  level  is  related  to  the  vertical  directional 
spectra  of  the  sea  surface -caused  sound. 
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1.0  INTRODUCTION 

The  experiment  described  in  this  technical  report  was  based  on  measure¬ 
ments  of  opportunity.  The  experiment  was  not  specifically  designed  to  meet 
the  objectives  stated  below.  Rather,  a  data  set  became  available  which 
appears  to  offer  additional  insight  into  the  space-time  characterization  of 
"sea  surface  sound"  generated  acoustic  ambient.  Such  a  statistical  character¬ 
ization  of  an  important  source  of  natural  underwater  acoustic  noise  has  clear 
relevance  to  the  design  of  systems  intended  to  detect  and  measure  underwater 
acoustic  radiation.  The  analysis  of  the  data  set  was  undertaken  to  meet  the 
objective  of  determining  the  wind  speed  dependence  of  high  frequency  vertical 
directional  spectra  caused  by  "sea  surface  sound." 

Significant  literature1"*  has  been  devoted  to  this  subject;  however, 
there  is  no  established  physical  theory  of  surface  noise  generation  at 
frequencies  above  approximately  50  Hz.5,3  "Bubbles  and  spray"  theories 
applicable  to  the  frequencies  of  interest  exist7"0  but  have  not  been  ade¬ 
quately  tested  or  integrated  into  space- time  statistical  descriptions.  Below 
the  wind  speed  for  which  incipient  whitecapping  exists,  there  is  even  less 
known.10  The  historical  approach  to  anisotropic  ambient  modeling  is  a  "source 
plane"  description11  derived  from  hydrophone  array  measurements.12  Thus  the 
modeling  efforts  are  data  driven  because  of  the  lack  of  a  basic  source 
description,  and  the  data  base  Is  flawed  because  much  of  the  data  was  taken  in 
open  ocean  areas  where  shipping13  and  sea  conditions  well  removed  from  the 
measurement1*  dominated  the  vertical  anisotiopic  field.  In  this  experiment 
external  contamination  of  the  data  is  controlled  by  the  location  of  the  exper¬ 
iment,  The  Tongue  of  the  Ocean  (TOTO)  in  The  Bahamas,  and  the  high  frequencies 
measured  (8  to  64  kHz).  The  location  is  a  relatively  deep  (1830  meters)  body 
of  water  isolated  from  global  shipping, 

The  next  section  of  this  report  gives  a  theoretical  view  of  the  role  of  a 
near-surface  layer  of  monopoles,  in  contrast  with  the  more  familiar  surface 
distribution  of  dipoles,  in  determining  the  acoustic  directional  spectra. 

That  section  is  followed  by  a  description  of  the  experiment  and  the  data 
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2 . 0  THEORY 

In  this  section  analytical  expressions  are  developed  for  the  vertical 
directional  spectra  caused  by  a  diffuse  layer  of  monopole  sources  beneath  an 
arbitrarily  rough  "pressure  release"  surface.  The  density  of  the  monopole 
sources  is  assumed  to  be  uniformly  distributed  in  the  horizontal  direction, 
and  either  exponential  or  uniform  distributions  are  considered  in  the  vertical 
direction.  The  objective  is  to  predict  the  change  in  the  horizontally 
directed  acoustic  field  due  to  a  change  in  the  near-surface  acoustic  source 
distribution.  The  effort  explores  the  hypothesis  that  the  origin  of  the 
acoustic  ambient  changes  from  momentum  fluctuations  on  the  boundary  when 
whitecaps  are  present  to  bubble  oscillation  beneath  the  surface  in  the  absence 
of  whitecaps. 

We  wish  to  model  a  monopole  source  at  a  depth  z0  beneath  the  surface  as 
pictured  in  figure  la;  i.e.,  the  positive  z  axis  is  below  the  surface.  At  a 
later  point  we  will  consider  a  monopole  density  function.  However,  here  we 
express  the  pressure  field  of  a  single  monopole  source  in  the  vicinity  of  the 
monopole,  i.e.,  at  a  radial  distance  R,  as 

P  “  P0/R  .  (1) 

where  P0  -  Q/4jt  and  Q  is  the  monopole  source  strength.  The  received  pressure 
at  a  large  distance  from  the  monopole  source  is 


P 

0  r 


ikR 


e  ikR' 
7  R,e 


-iwt 

e 

0  0 


(2) 


where  the  real  variable  n  is  the  water  attenuation  coefficient,  7  is  a 
complex-valued  reflect.'1  on  loss  coefficient,  and  k  is  the  wavenumber  magnitude 
in  the  radial  direction  of  propagation.  It  is  useful  to  write  the  surface 
reflection  coefficient  in  polar  form  as 
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7  “  Me 


-i2kf 


COB<f> 


(3) 


where  /j  is  an  amplitude  reflection  coefficient15  and  (2kf  cos <f>)  is  a  phase 
change18  from  a  reflection  from  a  "rough"  surface  of  height  relative  to  the 
mean  height  of  f  -  0. 


The  mean- square  received  pressure  becomes 
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where  the  bar  denotes  time  average,  and  AR  ~  R  -  R'  , 

Note  that  equation  (4)  will  be  realized  only  if  the  receiver's  spatial 
resolution  is  unable  to  separate  the  two  paths.  If  the  receiver  is  capable  of 
resolving  the  two  paths,  then  the  strengths  of  the  monopole  and  its  image  will 
be  observed  separately.  The  separability  of  the  two  arrivals  is  a  function  of 
the  ratio  of  the  source  and  receiver  depths  (z0/h) .  For  the  case  of  interest, 
for  which  the  receiver  depth  is  much  greater  than  the  maximum  source  depth 
(z0/h  «  1),  the  receiver  observes  the  monopole  and  image  arrival  from  the 
sauie  angle  but  with  a  phase  delay  (kAR)  .  Figure  1  shows  that  if  <j>  “  </>'  then 

fh 

^  cos^(z)  ’  anc*  (3a) 
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where  we  have  neglected  refraction  effects  only  over  the  short  spatial  inter¬ 
val  -z0  to  z0.  Thus  while  R'  ■  R  we  have 

AR  -  2z0  cos<£,  .  (6) 

In  light  of  these  approximations,  equation  (4)  becomes 


The  amplitude  reflection  coefficient15  is 


(7) 


■  2  2  -2k O  cold 

7  “  /i  -  e 


(8) 


where  a  is  the  root-mean-square  surface  roughness.  We  assert  that  the  random 

phase  and  amplitude  components  of  the  complex  reflection  coefficient,  equation  # 

(4),  are  statistically  independent,  and 


7  -  pe 
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for  a  Gaussian  distributed  surface  elevation.16  Thus,  by  using  equation  (9), 
equation  (7)  becomes 
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For  a  perfectly  smooth  surface  this  becomes  the  familiar  "dipole" 
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If  there  are  NM(z)  statistically  independent  monopoles  per  unit  volume,  then 


P2  -  Nm(z)P2 


(12) 


2 

is  the  monopole  source  pressure  level  per  unit  volume  and  dP£(z)  is  the  re¬ 
ceived  mean- square  pressure  level  per  unit  source  volume  at  range  R  from  the 
receiver. 


The  differential  received  mean-square  pressure  level  due  to  a  differen 
tial  source  volume  (r  dr  d 6  dz)  located  at  r,  6,  and  z  is 


dpz(z)  -  pz  r  dr  d$  dz 

r  r 

and  the  receiver  directional  spectral  function  is 


(13) 


p  r  dr  d(  —  , 
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dip  r  dip 


(14) 
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where  c is  the  differential  solid  angle  measured  at  the  receiver,  and  the 
differential  volume  (dA  dz)  is  measured  at  the  source.  This  problem  histori¬ 
cally  has  been  solved15  for  a  source  "plane"  which  involves  only  the  differen¬ 
tial  area.  Each  differential  source  area  contributes  to  a  unique  receiver 
elevation  angle.  For  the  "volume"  source  case  there  is  no  longer  a  unique 
relation  between  each  differential  source  volume  and  a  receiver  elevation 
angle.  This  point  is  most  easily  illustrated  for  the  isovelocity  case  for 
which  4>  -  <f>a  and 

tan^  -  — - -  for  tt/2  sS  <f>  <  ?r  . 

h  -  z 

o 

Thus  for  a  "source  plane"  both  h  and  z0  are  constant  and  r  uniquely  deter¬ 
mines  4>,  In  general  there  is  a  differential  "tube"  of  sources  which  contrib¬ 
ute  to  a  specific  elevation  angle.  Tha  total  acoustic  intensity  associated 
with  a  specific  direction,  9  and  <f> ,  is  determined  by  integrating  along  the 
tube,  i.e.,  incoherently  summing  all  of  independent  acoustic  sources  contrib¬ 
uting  to  a  specific  direction.  Within  the  limits  of  our  deep  receiver -shallow 
source  assumption,  i.e.,  D/h  «  1,  all  of  the  sources  in  the  volume 


fD 

dz  r  dr  d 9 

>  0 


contribute  to  a  unique  direction,  Thus  we  may  express  the  desired  directional 
spectra  as 


N M<*,*) 


rD 


dz 


J0 


dPr(z) 

6i> 


fD 

dz  Pr(z)  dfl  r  dr 
d<^  69  sin <j> 

J0 


and 


N  (  6  ,<t>) 

n 


rD  _ 

Pr(r,  0,  z)  dz 


“  0 


sind  d 6 
r  dr 


(15a) 


(15b) 
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Again  note  that  6  and  4>  are  measured  at  the  receiver.  Within  the  limits  of 
our  D/h  «  1  assumption,  r,  6,  and  R  arc  independent  of  z  and  thus 


N  <M> 

M 


A 6 

rD 

dz  N  (z) 

M 

1  +  a2  -  2a  cos  2kz  cos<t>  I 

V  J  V 

UJ  J 

0 

sind  dd 
r  dr 


(16) 


where 


a  -  e 


2  2 
-k  a 


00  4> 


The  source  and  receiver  angles  are  related  by  Snell's  law,  cos^I/c(z0)  - 
cos«^/c(h)  ,  where  c  is  the  fluid  phase  velocity.  Several  studies17,18  have 
seen  an  exponential  distribution  of  bubble  density  with  depth.  Thus  we 
evaluate  the  integral  in  equation  (16)  for  N(z)  -  N0e'l/D  after  extending  the 
integration  limit  to  +  ».  Thus 


N m(M) 


2 

■ 

r  ^ 

[P 

-2nR 

i  .  zl 

0 

u, 

o 

li+“j 

sln^  d^ 
r  dr 


2a 


f2kD  cos <t> 


(17) 


It  is  also  instructive  to  evaluate  equation  (16)  for  a  uniform  distribution  of 
monopoles,  N(z)  -  N0  for  z  £  D  and  zero  elsewhere.  In  this  case  we  get 


(p 

-2nR 

e 

(  2) 

sinf2kD  cos^  ' 

ll  +  a2 

.  2a - - V- 

UJ 

v.  ) 

^2kD  cos^  j 

sirnA  d  6 
r  dr 


(18) 
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For  Che  case  of  a  -  1,  equation  (18)  reduces  to  the  "volumetric  source" 
expression  of  Bannister,  Burgess,  and  Kewley.lfl 


There  are  numerous  physical  processes  at  the  sea  surface  which  cause 
momentum  fluctuations20  that  are  accurately  modeled  as  surface  dipole  sources. 
Analogous  with  equation  (16)  we  have 


N  d(d,<j>) 


2 

e"2nRN 

d 

2  , 

COS  <p 

a 

sind 

r 

M 

dr 

(19) 


where  NdPd2  is  the  dipole  source  strength  per  unit  surface  area. 


Before  evaluating  equations  (17),  (18),  and  (19) 
velocity-depth  profiles,  it  is  useful  to  evaluate  the 
isovelocity  ocean,  i.e,,  c(z)  -  c  and  n  -  0,  In  this 


for  arbitrary  sound 
equations  for  a  lossless 
case 


R2  -  r2  +  h2  -  r2  (l  +  cot2J  -  (^]  , 

<j>  -  <f>  ,  and 

S 

dd  cos 26  , 
dr  h 


and  for  surface  dipoles 


N  (6  ,<f>)  -  N  P2  cos 4> 
d  d  d 


(20) 


For  an  exponential  distribution  of  subsurface  monopoles  we  have 


N  (6  ,4>)  -  N  P2D 
M  0  0 


_ Jla _ 

1  +  (2kD  costf) ) 2 


COS  <t> 


(21) 


10 


TR  8537 


and  for  a  uniform  distribution  of  subsurface  monopoles  we  have 


nm(M)  -  nop;d 


(l  +  Q2) 


2a 


sin(2kD  c oaS) 


(2kD  coa<j>) 


cos  6 


(22) 


As  <f>- tt/2,  both  equations  (21)  and  (22)  become  dipole  sources  with  a 
frequency  (wavenumber)  dependent  strength 


lim 


N  (S  ,<t>)  -  2N  P2D(2kD)2 
M  0  0 


COSi f> 


(23) 


# 


This  result  is  independent  of  cr. 

It  is  useful  to  note  that  for  o  equal  to  zero 

lim  N  (8  ,<(>)  8N  P2R2k2D3  cos^  (dipole)  (24a) 

D— 0  M  000 

, .  2N  P2D 

^  V9"»  -  W  (monopole)  (24b) 

Thus,  as  the  depth  of  the  layer  of  monopoles  becomes  small,  the  direc¬ 
tional  spectra  becomes  the  same  as  that  due  to  a  "plane"  of  independent 
dipoles  with  a  vanishingly  small  intensity.  On  the  other  hand,  if  the  layer 
becomes  deep,  then  the  directional  spectra  becomes  the  same  as  that  due  to  a 
"plane"  of  independent  monopoles  whose  intensity  increases  with  the  first 
power  of  the  layer  depth. 

Thus  the  shape  of  the  vertical  directional  spectra  changes  significantly 
with  the  depth  of  the  layer.  From  equation  (16)  we  see  that  the  directional 
spectra  is  a  function  of  the  integral  of  the  monopole  depth  distribution. 

Next  we  are  interested  in  evaluating  the  directional  spectra  for  an 
exponentially  distributed  surface  layer  of  monopoles,  equation  (17),  and  a 
sound  velocity-depth  profile  typical  of  the  experiment  conditions  considered 


t 
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later  in  this  report,  Over  the  duration  of  the  experiment,  which  exceeded  1 
year,  a  wide  range  of  sound  velocity -depth  profiles  existed;  however,  because 
the  125 -meter  receiver  depth  nearly  always  exceeded  the  seasonal  thermocline 
depth,  the  source-receiver  "cross-duct"  configuration  was  typical.  In  the 
following  calculations  we  use  the  sound  velocity-depth  profile,  shown  in  fig¬ 
ure  2,  which  has  a  43-m  seasonal  thermocline  depth.  This  is  not  a  typical 
profile  but  rather  is  a  summer  profile.  It  was  chosen  because  it  reduced  the 
computation  load  of  the  calculations  because  the  amount  of  surface  contribut¬ 
ing  to  the  receiver  was  limited  by  refraction.  The.  "limiting  ray"  at  the 
receiver  due  to  a  near-surface  source  is  approximately  95°.  Note  that  in  the 
coordinate  system  of  figure  1,  which  has  the  positive  z  axis  pointing  down,  0 
is  down  and  n  is  up.  Because  of  the  hig,  r  quency  of  the  experiment  (8  to  64 
kHz)  we  use  eikonal  equations  to  model  the  -levant  acoustic  propagation  and 
we  treat  the  water  column  as  infinitely  deep.  This  assumption  will  be  justi¬ 
fied  in  the  next  section.  Figure  3  illustrates  the  dependence  of  the  direc¬ 
tional  spectra  on  the  depth  scale  of  the  exponential  layer  nondimens ionalized 
by  acoustic  wavelength.  The  calculation  is  performed  at  10  kHz  (15.5-cm  wave¬ 
length)  .  The  figure  illustrates  that  a  thin  layer  of  monopole  sources  having 
an  exponential  length  scale  exceeding  1  to  2  cm  will  have  a  10  kHz  directional 
spectra  observably  different  from  a  surface  dipole  spectra.  As  the  layer 
thickness  increases,  the  directional  spectra  moves  from  a  "dipole-like"  to  a 
"monopole-like"  pattern.  This  result  is  a  consequence  of  the  directional 
spectra  being  due  to  the  independent  addition  of  dipoles  (monopole  plus  nega¬ 
tive  image)  having  multiple  lobes.  The  deeper  the  layer,  the  more  multiple- 
lobe  contributors  present,  as  described  by  Hallett  and  Heitmeyer.21  From 
figure  3  we  conclude  that  once  the  exponential  length  scale  exceeds  1  wave¬ 
length,  the  directional  spectra  is  best  modeled  as  a  surface  layer  of  indepen¬ 
dent  monopole  sources.  Experimentally  then,  the  near-horizontally  directed 
acoustic  ambient  should  be  a  sensitive  observation  of  whether  monopole- like  or 
dipole-like  behavior  exists. 
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Figure  2.  Summer  Sound  Velocity-Depth  Profile 
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Figure  3.  Vertical  Directivity  Function  for  an  Exponentially  Distributed 

Monopole  Subsurface  Layer 
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Figure  3  illustrates  the  monopole -like  behavior  of  the  dipole  sources  as 
41  the  separation  of  the  elements  of  the  source  and  its  image  becomes  a  signifi¬ 

cant  fraction  of  a  wavelength  even  for  a  perfectly  reflecting  and  flat  sur¬ 
face,  i.e.,  an  unmodified  image  source.  Figure  4  illustrates  that  surface 
roughness  also  amplifies  the  monopole- like  behavior.  Physically  this  results 
41  from  a  modification  of  the  image  source,  as  would  be  expected.  Figure  5  pre¬ 

dicts  that  with  the  sound  velocity-depth  profile  of  figure  2  the  media  volume 
attenuation  will  also  modify  the  directional  spectra  at  16  kHz  and  above. 
Refractive  effects  limit  the  range  of  the  surface  contributing  to  the  ambient 
H  acoustic  field  to  less  than  about  2000  meters.  Above  16  kHz  the  attenuation 

loss  beyond  that  range  becomes  observable. 
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Figure  4.  Surface  Roughness  Effects  on  the  Vertical  Directivity  Function 
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3.0  EXPERIMENT  DESCRIPTION 


3.1  ACOUSTIC  DESCRIPTION 


Four  acoustic  sensors  were  used  in  this  experiment.  Three  of  them  were 
vertical  line  arrays  of  hydrophones  and  the  fourth  was  a  single  hydrophone. 
While  the  latter  was  certainly  the  least  "directive,"  it  too  had  some  direc¬ 
tional  characteristics  at  the  higher  frequencies  measured.  The  three  direc¬ 
tional  sensors,  i.e.,  line  arrays,  were  internally  wired  so  that  only  a  single 
output  was  available  from  each.  This  output  was  the  sum  of  the  individual 
hydrophones  of  the  array,  referred  to  as  the  "sum  beam."  The  line  arrays  were 
approximately  10  wavelengths  long  at  16,  32,  and  64  kHz.  The  data  in  this 
report  were  taken  at  8,  16,  32,  and  64  kHz.  The  8  and  16  kHz  data  were  taken 
from  the  array  designated  L1(  and  the  32  and  64  kHz  data  were  taken  from 
arrays  L2  and  L3.  The  calculated  vertical  pattern  function  for  each  of  these 
arrays,  which  contained  nonuniform  weighting,  are  pictured  in  figures  6 
through  9.  The  half-power  width  of  the  vertical  angular  window  was  approxi¬ 
mately  5°  at  16,  32,  and  64  kHz,  and  10°  at  8  kHz.  The  pattern  function  of 
each  hydrophone  in  the  arrays  was  azimuthally  uniform  and  mildly  nonuniform  in 
elevation  angle.  The  single  hydrophone  was  a  composite  cylinder  having  a  ver¬ 
tical  axis.  As  was  the  case  with  the  arrays,  the  hydrophone  had  an  azimuth- 
ally  uniform  response  at  all  frequencies.  In  the  vertical  direction  the 
hydrophone  pattern  function  was  mildly  directional  at  16  kHz  and  significantly 
so  at  32  and  64  kHz. 

Measured  patterns  exist  for  the  hydrophone.  The  patterns  are  rather  com¬ 
plex  in  the  "low  response"  directions  but  fairly  uniform  in  the  "high 
response"  directions.  It  was  useful  to  model  the  hydrophone  response  in  such 
a  way  that  the  model  fit  the  "high  response"  directions  well  and  fit  the  "low 
response"  directions  less  well.  The  requirement  of  the  model  was  to  ade¬ 
quately  predict  the  total  hydrophone  output  in  measuring  acoustic  ambient 
fields.  The  model  was  created  by  taking  the  product  of  a  "ring"  transducer 
and  a  "line"  transducer  modeled  as  a  zero-order  Bessel  function  and  a  "sync" 
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Figure  6.  Array  Lx  Vertical  Pattern  Function  at  8  kHz 


Figure  7.  Array  Lj  Vertical  Pattern  Function  at  16  kHz 
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function,  The  coefficients,  i,e. ,  nondimens ional  cylinder  radius  and  length, 
were  modified  to  "best"  fit  the  measured  patterns,  The  explicit  algorithm  is 


where  J0  is  a  zero-order  Bessel  function,  <£  is  the  elevation  angle,  X  is  the 
acoustic  wavelength,  and  d  and  L  are  the  modified  cylinder  diameter  and  length 
scales,  Figures  10  through  13  show  the  calculated  hydrophone  vertical  pattern 
function  at  the  four  frequencies  of  interest. 
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Figure  10.  Omnihydrophone  Vertical  Pattern  Function  Model  at  8  kHz 


Figure  11.  Omnihydrophone  Vertical  Pattern  Function  Model  at  16  kHz 
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Figure  12.  Omnihydrophone  Vertical  Pattern  Function  Model  at  32  kHz 
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Figure  13.  Omnihydrc phone  Vertical  Pattern  Function  Model  at  64  kHz 


22 


3.2  EXPERIMENT  GEOMETRY 


The  acoustic  sensors  were  located  in  The  Tongue  of  the  Ocean  (TOTO) , 

The  Bahamas.  They  were  positioned  in  the  middle  of  the  major  axis  of  the 
tongue  at  a  latitude  of  24°40'  N.  The  sensors  were  at  a  depth  of  125  m  in 
water  1830  m  deep.  The  location  was  approximately  12  km  from  the  nearest 
escarpment.  If  we  define  a  "range  curtain"  to  be  the  range  at  which  the 
geometric  spreading  loss  is  matched  by  the  attenuation  loss,  then  the  maximum 
range  of  influence  was  less  than  35  km  for  frequencies  greater  than  8  kHz. 

All  data  were  examined  for  signs  of  extraneous  acoustic  signals,  and  all.  sus¬ 
pect  data  sets  were  discarded. 

The  four  sensors  were  moored  to  the  bottom,  and  each  of  the  sensors  was 
electrically  hard  wired  to  the  AUTEC,  Andros  Island  Site  1  facility  for 
recording  and  analysis.  Data  were  available  from  20  May  1988  to  3  April  1989. 
During  that  period  215  data  sets  of  appropriate  quality  became  available  for 
this  analysis.  The  approximate  1-year  period  of  the  experiment  sampled  all 
seasons  of  the  local  sound  velocity-depth  profile  shown  in  figure  14.  The 
125 -m  depth  of  the  sensors  means  that,  aside  from  the  short  period  for  which 
the  winter  profile  existed,  the  measurements  were  made  beneath  the  seasonal 
thermocline . 

3.3  DATA  ANALYSIS 


The  analysis  described  in  section  2  treats  a  "bottomless"  ocean,  i.e., 
only  down-going  acoustic  propagation  is  possible.  To  check  the  conformity  of 
this  approach  with  the  realities  of  the  experiment,  the  directional  spectra 
were  estimated  using  the  techniques  described  by  Kennedy  and  Szlyk  in  1988. 22 
The  approach  features  a  propagation  loss  algorithm  that  obtains  full  solutions 
of  the  acoustic  wave  equation,  surface  roughness,  arbitrary  sound  velocity- 
depth  profiles,  arbitrary  bottom  geoacoustic  parameterizations ,  and  a  plane  of 
independent  dipole  radiators  located  on  the  water  surface.  An  important  weak¬ 
ness  of  the  approach  is  that  it  ignores  all  range  dependency.  In  the  geometry 
of  this  experiment,  range -dependent  bathymetry  is  an  issue.  Thus  the  predic¬ 
tion  ignores  steep- to- shallow  angle  direction  conversion  that  might  occur  at 
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Figure  14.  Sound  Velocity-Depth  Profile  of  TOTO 


the  escarpments.  The  results  of  the  calculations  are  shown  in  figures  15  and 
16  which  illustrate  the  predictions  for  summer  and  winter  profiles,  i.e., 
characterized  by  the  depth  of  the  seasonal  thermocline  (figure  14).  The  spe¬ 
cific  bottom  loss  characterization  used  for  the  8  through  64  kHz  frequency 
range  is  from  Marsh  and  Schulkin.23  The  figures  indicate  the  reflected  bottom 
energy  in  this  frequency  range  is  in  excess  of  15  dB  less  than  the  direct 
energy.  The  conclusion  is  that  a  bottomless  ocean  is  a  reasonable  model. 

The  autospectral  density  of  each  of  the  four  sensors  was  determined  using 
constant-bandwidth  Fourier  transform  analysis  techniques.  The  spectral  reso¬ 
lution  of  the  measurements  using  a  Hanning  window  was  15,  37.5,  and  75  Hz  for 
the  three  analysis  ranges.  Fifty  nonoverlapped  transforms  were  averaged  to 
estimate  the  spectral  density.  The  confidence  interval  of  these  spectral 
estimates  is  +1.1  dB.  The  spectral  estimates  were  compared  with  anticipated 
electronic  noise  limitations.  This  was  not  a  straightforward  operation 
because  full  system  electronic  noise  "floors"  were  not  readily  obtained;  how¬ 
ever,  any  data  failing  to  meet  the  anticipated  voltage  limitations  were  dis¬ 
carded.  The  autospectral  density  values  from  each  of  the  sensors  at  8,  16, 

32,  and  64  kHz  were  then  manually  recorded. 

The  line  antenna  array  gain  was  evaluated  by  taking  the  ratio  of  the  line 
array  output  to  the  hydrophone  output.  The  array  gain  was  then  compared  with 
a  calculated  directivity  index,  to  be  described  shortly,  and  the  difference, 
which  was  labeled  anisotropic  gain,  was  recorded  in  a  data  base  that  contained 
wind  speed  measurements  made  at  the  time  of  the  array  output  measurement.  The 
hydrophone  autospectral  density  was  also  included  in  the  data  bjse.  The  ane¬ 
mometer  used  was  located  10  m  above  the  sea  surface  on  a  small  tower  located 
in  the  shallows  adjacent  to  Andros  Island.  The  anemometer  was  approximately 
13  km  from  the  acoustic  sensor  location.  The  data  base  was  accumulated  on  a 
mainframe  computer  which  then  was  used  to  generate  the  results  described  in 
the  next  section  of  this  report, 
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Frequency-Directivity  Spectra  for  a  Summer  Dipole  Noise  Model 
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Frequency-Directivity  Spectra  for  a  Winter  Dipole  Noise  Model 
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The  directivity  index  that  served  as  the  reference  for  the  hydrophone  line 
array  outputs  was  calculated  so  that  it  represented  the  directivity  index  of 
the  measurement  hardware.  Thus  all  the  acoustic  antenna  results  were  normal¬ 
ized  by  the  directivity  index,  which  was  found  by  calculating  the  ratio  of  the 
outputs  of  each  of  the  three  arrays  to  the  output  of  the  reference  hydrophone, 
all  in  the  presence  of  a  three-dimensional  isotropic  pressure  field.  The  cal¬ 
culation  was  made  using  the  equation 

B(<j&,f)H  ,  f )  sin^  d<t> 

^ _ - _  ,  (26) 

|  H  (tf> ,  f)  sin^  &<t> 

J0  2 

where  B(^,f)  is  the  pattern  function  of  the  appropriately  weighted  and  spaced 
"point"  hydrophones,  H1(^,f)  is  the  pattern  function  of  the  individual 
hydrophones  used  in  the  antennas,  and  H2(^,f)  is  the  directional  pattern  of 
the  hydrophone  calculated  in  equation  (25)  .  Correcting  measured  data  to 
account  for  sensor  pattern  functions  is  traditionally  difficult  because  the 
data  represent  the  integral  of  the  product  of  the  pattern  function  and  the 
solid  angle  density  function  of  the  field  variable  being  measured.  Thus  the 
experimenter  is  in  the  dilemma  of  needing  the  result  of  the  measurement  prior 
to  making  the  correction.  This  was  avoided  here  by  not  "correcting"  the  array 
gain  measurement  but  rather  by  comparing  it  with  the  directivity  index  which, 
by  definition,  assumes  an  explicit  solid-angle  density  function,  However,  in 
presenting  the  hydrophone  output  as  a  measure  of  the  total  acoustic  intensity, 
it  is  necessary  to  make  a  correction  for  the  transducer  vertical  pattern  func¬ 
tion.  We  will  make  the  correction  by  letting  Sh(w)  be  hydrophone  (measured) 
autospectral  density 

S .  («>  -  [  dtf  H  <w,0)N(«,tf)  ,  (27) 

h  J4w  2 

where  \p  is  the  solid  angle,  H2  is  the  hydrophone  pattern  function,  and  N  is 
the  actual  noise  field  directional  function.  An  omnidirectional  hydrophone 
would  measure 
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Forming  the  ratio 


f  dip  H  (w,ip)N(u,ip) 
(a ))  J4m  Z 


,<“>  f  c 

J4jt 


(29) 


d^  N(w,i/>) 


requires  an  explicit  form  for  N(u,^i).  Following  Cox,24  we  expand  the 
directional  spectra  in  terms  of  Legendre  polynomials,  Pn(cos<6),  after 
asserting  azimuthal  symmetry  of  the  pressure  field.  The  result  is 


N(w,^)  -  X  cn(w)Pn(cos^)  . 
n-0 


(30) 


The  advantage  of  this  form  is  that  the  field  is  expressed  as  perturba¬ 
tions  (n  >  0)  on  an  isotropic  field  (n  -  0) .  This  follows  from 

Po(cos0)  -  1  . 


(31) 


In  correcting  for  the  hydrophone  pattern  function  in  our  measurements  we 
limit  ourselves  to  a  first-order  isotropic  approximation  so  that 


S  (w)  —  S  («) 
o  h 


f 


sin<£  d<p 


(32) 


f 

JQ 


sin^  dip  H2(w,<^) 


Attention  needs  to  be  drawn  to  the  ''hole1'  in  the  directivity  spectra 
(at  ip  -  ?r/2)  that  is  predicted  when  the  receiving  hydrophones  are  below  the 
seasonal  thermocline.  This  "hole"  is  seen  in  figure  15,  as  well  as  in  figures 
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3  through  5.  The  sound  velocity-depth  profile  used  to  generate  these  figures 
used  a  fast  surface  velocity  (1545  m/s)  which  maximized  the  size  of  the 
"hole."  In  general  the  size  of  "hole,"  relative  to  the  size  of  the  half -power 
beamwidth,  is  frequency  and  sound  velocity  dependent,  so  the  frequency  varia¬ 
tion  of  a  beamformer  output  would  be  expected  to  decrease  as  the  frequency 
increases  to  the  point  that  the  beampattern  half-power  beamwidth  becomes  nar¬ 
rower  than  the  "hole."  The  beamformer  output  usually  decreased  at  approxi¬ 
mately  the  3  dB  per  octave  rate  one  would  expect  without  the  "hole"  present. 
Thus  the  "hole"  was  not  identified  in  the  data. 

The  two  fundamentally  independent  variables  of  the  experiment  are  wind 
speed  and  frequency.  The  measured  wind  speed  at  a  10-m  height  is  better 
related  to  "sea  surface  sound"  by  converting  the  measurement  to  friction  velo¬ 
city  (U*) ,  i.e.,  atmospheric  stress  on  the  ocean  surface,  as  done  by  Kerman.7 
The  published  algorithm  of  Wu2S  was  used  to  make  the  conversion.  Friction 
velocity  will  be  presented  in  a  nondimens ional  form  by  normalizing  the  minimum 
phase  velocity  (U0)  of  the  capillary- gravity  wave  field.7  A  unity  value  of 
the  nondimens ional  friction  velocity  will  be  used  as  the  point  of  incipient 
whitecapping. 
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4.0  EXPERIMENT  RESULTS 


4.1  ACOUSTIC  INTENSITY 


# 


The  hydrophone  autospectral  densities  are  displayed  as  a  function  of  wind 
speed  in  figures  17  through  20.  These  data  are  compared  with  previous 
measurements  to  establish  a  historical  relation.  All  data  sets  are  compared 
with  the  well-known  globally  averaged  Knudsen  curves.2  In  addition,  the  8  and 
16  kHz  data  could  be  compared  with  the  Caron  and  Sheffield  data.26  The  latter 
data  were  taken  in  the  same  area  some  years  earlier.  One  sees  that  the  pre¬ 
sent  data  are  consistent  with  the  Caron  and  Sheffield  data,  as  they  should  be, 
and  are  similar  to  but  consistently  lower  than  the  Knudsen  data  except  at  64 
kHz.  The  lack  of  predicted  wind  speed  dependence  in  the  64  kHz  data  is  unex¬ 
pected,  and  it  can  only  be  assumed  that  the  sensor  was  not  operating  correctly 
at  that  frequency,  possibly  because  of  an  undetected  electronic  noise  limita¬ 
tion.  No  further  analysis  of  the  64  kHz  data  was  undertaken. 

Figure  21  illustrates  the  frequency  dependence  of  the  data.  The  figure 
was  determined  by  grouping  the  data  into  the  wind  speed  cells  shown  in  the 
figure  at  each  of  the  selected  frequencies.  A  regression  line  was  than  fit 
across  the  frequency  to  each  wind  speed  cell.  The  resulting  "power  law"  fit 
to  the  data  indicated  a  frequency  dependence  that  varied  from  -1.74  at  the 
lowest  wind  speed  cell  to  -2.1  at  the  highest  wind  speed  cell.  The  slopes  are 
comparable  with  the  "classical"  5  to  6  dB/octave,  i.e.,  -1.7  to  -2.0,  associ¬ 

ated  with  the  Knudsen  curves.  These  same  slopes  are  obtained  from  physical 
arguments.  Measured  acoustic  spectra27  due  to  spray  on  the  surface  is  -1.7. 

A  -2  dependence  was  obtained  by  Prosperetti  on  purely  theoretical  grounds  when 
modeling  the  acoustic  radiation  of  oscillating  bubbles  resulting  from  a  white- 
cap  . 28 

The  previous  hydrophone  sound  pressure  level  data  were  replotted  as  a 
function  of  nondimensional  friction  velocity,  i.e.,  surface  friction  velocity 
nondimensionalizad  by  the  minimum  phase  speed  of  the  capillary-gravity  wave. 
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Figure  21. 


Frequency  Dependence  of  SPL  Data  at  Various  Wind  Speeds 
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The  traditional  idea  is  that  a  value  of  1  for  this  nondimens ional  variable 
represents  incipient  whitecapping.  Each  of  the  data  sets  is  partitioned  into 
"whitecaps  present"  (WP)  and  "whitecaps  not  present"  (WNP)  according  to 
whether  the  nondimens ional  friction  velocity  is  greater  than  or  less  than  1, 
respectively.  A  regression  line  is  fit  to  each  partition.  The  results  for 
each  frequency  are  presented  in  figures  22  through  24.  The  slope  of  the  line 
is  seen  to  be  different  for  the  WP  and  WNP  cases.  For  the  WP  case  the  power 

law  exponent  varies  from  2.2  to  2.9.  For  the  WNP  case  the  exponent  varies 

from  1.2  to  1.9.  A  composite  of  all  of  the  data  is  obtained  by  adjusting  the 
curves  so  that  the  center  of  mass  of  the  data  at  a  nondimens ional  friction 
velocity  of  1  coincides.  This  is  illustrated  in  figure  25.  The  resulting 
power  law  fit  to  the  composite  data  yields  2.5  and  1.6  respectively  for  the 
WNP  and  WP  cases,  This  is  a  somewhat  smaller  exponent  for  the  WNP  case  than 

the  3.0  reported  by  Kerman.7  However,  for  the  WP  case  the  agreement  is  quite 

good,  i.e.,  1.5  vice  1.6.  There  appears  to  be  no  physical  argument  for  either 
power  law.  One  obvious  characteristic  of  the  data  being  presented  is  the 
larger  variability  of  the  WNP  data  relative  to  the  WP  data.  This  observation 
is  quantified  in  figures  26  and  27  which  are  histograms  of  the  residuals 
plotted  to  the  regression  lines  of  figure  25  for  the  WNP  and  WP  cases.  The 
difference  is  striking. 

4.2  ANISOTROPIC  ARRAY  GAIN 


# 


We  explicitly  defined  a  term  in  section  3.3,  i.e.,  anisotropic  array 
gain,  to  describe  the  difference  in  the  realized  gain  of  the  array  relative 
to  the  directivity  index  of  the  array.  By  array  gain  we  mean  the  improve-  • 

ment  in  signal- to-noise  ratio  (SNR)  between  the  SNR  of  the  homogeneous 
pressure  field  measured  by  the  array's  hydrophones  and  the  array  beamformer 
output  when  "steered"  to  the  horizontal  direction.  Array  gain  is  measured 
here  without  a  signal  field  present  by  "calibrating"  the  beamformer  output  so  ^ 

that  its  response  tc  an  acoustic  plane  wave  arriving  from  the  "steered"  direc¬ 
tion  is  at  the  same  acoustic  level  as  seen  by  the  hydrophone.  In  this  case, 
which  is  typical  of  the  industry,  the  signal  gain  is  unity  and  the  measured 

• 
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NONDIMENSIONAL  FRICTION  VELOCITY 

Figure  24.  Omnihydrophone  SPL  Versus  Nondimensional 
Friction  Velocity  at  32  kHz 
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Figure  25.  Composite  SPL  Versus  Nondimensional  Friction  Velocity 
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noise  reduction  is  exactly  the  reciprocal  of  the  array  gain.  We  then  calcu¬ 
lated  a  term  which  we  call  the  anisotropic  array  gain  which  measures  the 

array  gain  due  to  the  presence  of  vertical  anisotropy  in  the  acoustic  ambient 
field  by  normalizing  the  measured  array  gain  by  a  calculated  directivity 
index  (equation  (26)). 

Figures  28  through  30  show  that  the  anisotropic  gain,  when  "steered"  to 
the  horizontal  direction,  is  significantly  wind  speed  dependent.  In  general 
the  gain  exceeds  that  of  an  isotropic  noise  field  at  high  wind  speeds  and  is 
less  than  the  isotropic  gain  at  low  wind  speeds.  The  frequency  dependence  of 
these  observations  is  shown  in  figure  31,  which  illustrates  a  mild  frequency 

variation.  The  functional  fit  to  the  data,  grouped  by  wind  speed,  is  a 

second-order  polynomial  regression.  The  16  kHz  peaking  appears  to  be  physi¬ 
cally  significant.  To  examine  this  wind  speed  dependence,  we  again  use  the 
nondimens ional  friction  velocity  as  the  independent  variable  in  figures  32 
through  34.  As  was  the  case  for  the  acoustic  intensity,  there  is  a  clear 
change  in  the  character  of  the  curve  for  unity  value  of  the  independent  vari¬ 
able,  indicating  a  physical  process  change  between  the  whitecaps  present  case 
and  the  whitecaps  not  present  case.  A  nondimens ional  anistropic  gain  was 
formed  by  equating  the  center  of  mass  of  the  data  for  each  frequency  at  a  non- 
dimensional  friction  velocity  of  1.  This  nondimens ional  data  composite  is 
shown  in  figure  35.  Each  data  set  was  again  subdivided  into  the  two  cases  and 
a  power  law  was  fit  to  each  subset.  In  the  WNP  case  the  individual  power  var¬ 
ied  from  1.6  to  2.1,  with  the  composite  value  being  1.9.  In  the  WP  case  the 
individual  powers  varied  from  0.04  to  0.12,  with  the  composite  being  0.1, 

Thus  the  array  gain  in  the  horizontal  direction  is  nearly  independent  of  wind 
speed  if  whitecaps  are  present  and  changes  at  some  19  dB  per  decade  if  white- 
caps  are  not  present. 

As  was  true  for  the  acoustic  intensity,  the  spread  of  the  residuals 
about  the  regression  line  fit  is  noticeably  different  between  the  WNP  and 
WP  cases.  This  is  shown  clearly  in  figures  36  and  37,  which  are  histograms  of 
the  residuals  in  each  case.  However,  what  is  unique  to  the  anisotropic  gain 
case  is  the  nonsymmetric  distribution  of  the  WP  case  histogram.  This  visually 
apparent  characteristic  of  the  data  was  checked  by  calculating  a  measure  of 
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Figure  30.  Anisotropic  Gain  Versus  Wind  Speed 

at  32  kHz 
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figure  36.  Histogram  of  Anisotropic  Gain  Residuals  With  Whitecaps  Not  Present 


Figure  37.  Histogram  of  Arisotropic  Gain  Residuals  With  Whitecaps  Present 
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skewness.20  This  metric  varied  from  0.2  in  the  WNP  case  to  0.8  in  the  WP 
case.  While  no  formal  statement  of  statistical  significance  was  worked  out 
for  this  difference,  it  is  clear  that  the  WP  case  appears  to  have  a  physical 
upper  limit  and  that  all  variation  is  constrained  by  this  limit.  This  is  not 
true  in  the  WNP  case  which  has  symmetrical  variation  about  the  power  law. 

Let  us  return  to  the  initial  observation  that  the  variability  of  the 
data  in  the  WNP  case  is  greater  than  in  the  WP  case.  One  obvious  possibility 
is  that  data  contamination  due  to  an  electrical  noise  floor  at  the  lower  spec¬ 
tral  levels  of  the  WNP  case  is  contributing  to  the  variability.  Another  pos¬ 
sibility  is  that  without  the  uniform  generation  of  spray  and  bubbles  by  the 
whltecaps,  the  sources  of  acoustic  ambient  are  more  irregular  and  dependent  on 
an  increased  number  of  independent  variables .  Examples  would  be  the  distribu¬ 
tion  and  excitation  of  bubbles  due  to  the  previous  history  of  the  sea  condi¬ 
tions  as  well  as  the  present  state  of  the  current,  near- surface  turbulence, 
and  internal  waves.  An  indication  that  the  latter  cause  might  be  the  case  is 
illustrated  In  figures  38  to  40,  which  are  plots  of  anistropic  gain  with  total 
acoustic  intensity  as  the  independent  variable.  The  idea  is  that  the  total 

acoustic  intensity  is  a  measure  of  the  acoustic  state  of  the  sea,  where,  for  a 

given  wind  speed  (without  whltecaps  present) ,  a  range  of  acoustic  states  is 
possible  and  it  is  that  range  of  conditions  which  caused  the  large  variability 

being  discussed.  The  reduced  variability  in  the  data  at  low  acoustic  inten¬ 

sity  seen  in  the  figures  shows  that  total  acoustic  intensity  is  a  better  pre¬ 
dictor  of  the  anisotropic  gain  at  low  wind  conditions  than  was  wind  speed- 
derived  quantities.  If  the  change  in  acoustic  state  were  caused  by  sources 
which  would  be  contaminants  to  the  data,  such  as  shipping  or  biological 
sources,  then  the  resulting  anisotropic  gains  would  have  been  controlled  by 
the  sidelobes  of  the  antenna,  causing  increases  in  the  gain,  and  the  variabil¬ 
ity  would  have  increased  rather  than  decreased,  as  was  the  case,  In  the  next 
section  we  will  argue  that  the  range  in  acoustic  states  of  the  sea  might  be 
caused  by  the  range  in  bubble  layer  depth.  We  note  that  some  sound  velocity- 
depth  profile  effects  could  be  due  to  the  refractive  "hole"  discussed  in  sec¬ 
tion  3.  This  would  add  variability  to  the  anisotropic  gain  measurements  both 
in  the  WP  and  WNP  cases. 
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SPL  SPECTRA  (dB// ft  Pa^/Hz) 

Figure  40.  Anisotropic  Gain  Versus  Acoustic 
SPL  Spectra  at  16  32  kHz 


5.0  CONCLUSIONS 


In  this  section  we  will  draw  some  tentative  conclusions  about  the  causa 
of  the  high  frequency  acoustic  ambient  directivity  spectra  by  bringing 
together  the  results  of  the  theoretical  discussions  of  section  2  and  the 
measured  results  of  section  4.  The  conclusions  are  necessarily  tentative  for 
two  reasons:  (1)  the  experiment  only  measured  the  horizontally  directed 
acoustic  ambient  and  is  an  incomplete  observation,  and  (2)  even  complete 
knowledge  of  the  sound  field  is  not  sufficient  to  determine  uniquely  the 
scarce  as  demonstrated  by  Dowling  and  Ffowcs  Williams.30 

We  begin  the  discussion  by  reviewing  the  four  basic  observations  of  the 
experiment: 

1.  The  anisotropic  gain  is  a.  function  of  wind  speed,  as  illustrated  in 
figure  35. 

2.  The  functional  dependence  of  the  anisotropic  gain  on  wind  speed  is 
dichotomized  according  to  whether  or  not  whitecaps  are  present,  as  illustrated 
In  figure  35. 

3.  The  probability  density  function  of  the  measured  anisotropic  gain  at 
a  given  wind  speed  is  skewed  to  lower  values  with  whitecaps  present,  indicat¬ 
ing  a  physical  limit  to  the  gain,  as  shown  in  figures  36  and  37. 

4.  The  variability  of  the  anisotropic  gain  without  whitecaps  present  is 
significantly  greater  than  with  whitecaps  present,  and  this  variability  is 
reduced  by  substituting  total  acoustic  intensity  for  wind  speed  as  the  inde¬ 
pendent  variable. 

The  derived  quantity  in  the  theoretical  discussion  of  section  2  is  the 
directivity  spectra,  while  the  measured  quantity  of  section  4  is  the  line 
array  beamformer  output  steered  to  the  horizontal  direction.  The  well-known 
relation  between  these  two  quantities,  for  an  azimuthally  uniform  field,  is 
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S*(f)  -  2*  B(^,f)H  U,f)N(tf,f)  sin*  d* 

b  “  0  1 


(33)  • 


Thus,  we  can  expect  the  horizontally  directed  beamformer  output  to  be  a 
sensitive  measure  of  the  vertical  directional  spectra  in  the  horizontal 
direction.  From  section  2  we  saw  that  the  near-horizontal  values  of  the 
directional  spectra  are  the  most  physically  revealing  aspect  of  the  func¬ 
tion.  This  is  true  to  the  extent  that  the  discussion  is  complete,  i.e.,  the 
acoustic  ambient  diffuse  sources  consist  of  surface  dipoles  and  below-surface 
monopoles  only.  These  two  source  types  lead  to  quite  different  source  effi¬ 
ciencies  at  shallow  "launch"  angles.  While  the  dipole  is  known  to  have  a  null 
in  the  horizontal  direction,  the  layer  of  monopoles  leads  to  a  directional 
spectra  that  is  determined  by  the  integral  of  the  vertical  density  distribu¬ 
tion  weighted  by  a  factor  that  describes  coherency  of  the  surface  reflection. 
This  form  is  shown  explicitly  in  equation  (16).  Thus  a  layer  of  monopoles 
beneath  the  ocean  surface  can  asymptotically  lead  to  either  surface  monopole - 
like  or  surface  dipole-like  behavior.  The  surface  monopole-like  behavior 
results  from  contributions  of  individual  monopoles  that  are  greater  than  a 
significant  fraction  of  a  wavelength  from  the  sea  surface,  i.e.,  acoustically 
deep.  As  the  depth  of  the  individual  contributors  increases,  the  pattern 
function  of  the  monopoly  source  and  its  surface  reflection  produces  a  multi¬ 
pole  function  which,  in  concert  with  other  such  sources  at  incrementally  dif¬ 
ferent  depths,  leads  to  a  composite  monopole  structure.  Any  reduction  in  the 
coherency  of  the  surface  reflection  caused  by  a  non- smooth  function  also  pro¬ 
duces  a  composite  monopole  structure.  The  surface  dipole-like  behavior 
results  from  the  opposite  conditions  wherein  all  the  subsurface  monopoles  are 
loss  than  a  small  fraction  of  an  acoustic  wavelength  from  the  surface  (acous¬ 
tically  shallow),  and  the  surface  produces  only  npecular  reflections.  Guo  has 
pointed  out  that  in  this  case  the  resulting  dipole  has  a  vanishingly  small 
source  level.®  Changes  in  the  behavior  of  the  directional  spectra  of  a 
below-surface  layer  of  monopoles  is  best  observed  at  near-horizontal  angles, 
as  illustrated  in  figures  3  through  5, 


« 
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In  the  absence  of  surface  dipole  sources,  the  subsurface  acoustically  deep 
layer  of  monopole  sources  would  lead  to  the  situation  in  which  the  acoustic 
ambient  is  dominated  by  horizontally  directed  energy.  This  r' suit  would 
become  more  pronounced  as  the  layer  depth  and/or  surface  roughness  increases. 
Because  the  acoustically  thin  layer  of  monopoles  leads  to  a  vanishingly  small 
dipole  source,  it  is  of  little  practical  interest..8,20  However,  at  the 
frequencies  of  interest  in  this  report,  surface  dipole  sources  would  result 
from  sprays  of  droplets  on  the  surface , 20,7,7  The  surface  dipole  sources  are 
known  to  lead  to  the  situation  in  which  the  acoustic  ambient  is  dominated  by 
vertically  directed  energy.  Thus,  when  evaluating  the  anisotropic  gain  for 
horizontal  steering,  the  surface  monopole  and  surface  dipole  cases  lead  to 
distinctly  different  values.  In  the  former  case  one  is  11  looking  away  from  the 
noise,"  i.e.,  high  array  gain,  and  in  the  latter  case  one  is  "looking  into  the 
noise,"  i.e.,  low  array  gain.  Furthermore,  we  note  that  the  surface  dipole 
case  leads  to  a  maximum  gain  situation;  i.e. ,  the  maximum  possible  gain  in  the 
horizontal  direction  occurs  when  the  acoustic  ambient  is  caused  by  a  surface 
layer  of  dipole  sources. 

The  above  models  can  be  shown  to  fit  qualitatively  our  experimental  obser¬ 
vations  as  follows: 

1.  A  wind  speed- dependent  array  gain  would  be  expected  if  the  wind  and 
the  resulting  sea  changed  the  source  structure  from  surface  dipoles  to  a  sub¬ 
surface  layer  of  acoustically  deep  mcnopoles ,  presumably  oscillating  bubbles. 
Furthermore,  the  wind,  as  well  as  other  factors,  would  be  expected  to  change 
the  vertical  distribution  of  tne  monopole  density,  and  such  changes  would 
change  >.he  array  gain. 

2.  With  whitecaps  present  the  acoustic  sources  would  be  expected  to  be 
dipole- like  as  a  result  of  a  ready  supply  of  spray  and  acoustically  thin 
layers  of  bubbles.  These  dipoles  would  minimize  the  horizontal  values  of  the 
ambient  directional  spectra  so  that  the  horizontally  steered  array  gain  would 
exceed  the  directivity  function.  In  fact  there  would  be  a  maximum  value 
determined  by  a  pure  surface  dipole  source.  Thus  the  largest  values  of  the 
horizontally  steered  array  gain  would  be  expected  to  occur  v.'ith  whitecaps 


49 


TR  8537 


present,  as  the  measurements  Indicated.  Furthermore,  because  an  array  gain 
maximum  exists,  the  statistical  distribution  of  values  would  be  skewed  to 
lower  values,  as  was  measured. 

3.  With  no  whitecaps  present  the  supply  of  spray  and  acoustically  thin 
layers  of  monopoles  would  be  significantly  diminished.  In  the  absence  of 
these  sources  the  role  of  an  acoustically  thick  layer  of  bubbles  would  become 
more  important.  The  presence  of  this  layer  of  sources  would  increase  the  hor¬ 
izontal  values  of  the  directivity  function  relative  to  the  vertical  values  and 
the  array  gain  would  be  expected  to  be  lower  than  was  the  case  when  whitecaps 
were  present,  as  occurred  in  the  data.  Because  the  wind  would  not  be  the  only 
cause  for  the  vertical  distribution  of  the  monopoles,  one  would  expect  a  sig¬ 
nificant  range  of  conditions  to  exist  at  a  particular  wind  speed.  This  var¬ 
iability  was  certainly  present  in  the  experiment  results, 

4.  If  the  above  reasoning  is  correct,  then  the  data  argue  that  in  gen¬ 
eral  the  vertical  depth  scale  of  the  monopole  layer  increases  with  decreasing 
wind  speed  although  with  significant  variation  in  the  value  for  a  given  wind 
speed.  The  depth  scale  variation,  or  more  exactly  the  array  gain  variation, 
is  directly  related  to  the  total  acoustic  intensity  of  the  acoustic  ambient, 
which  serves  as  a  good  predictor  of  the  directional  spectra. 
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